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Measurement of the charged current quasi-elastic cross-section 
on carbon with the T2K on-axis neutrino beam 
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We report a measurement of the charged current quasi-elastic cross-sections on carbon in 
the T2K on-axis neutrino beam. The measured charged current quasi-elastic cross-sections on 
carbon at mean neutrino energies of 1.94 GeV and 0.93 GeV are (11.95 ± 0.19(stat.) t\-f^{syst.)) X 



10“®®cm^/neutron, and (10.64 ± 0.37(stat.)l^'gg(st/st.)) x 10“®®cm^/neutron, respectively. These 
results agree well with the predictions of neutrino interaction models. In addition, we investigated 
the effects of the nuclear model and the multi-nucleon interaction. 
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I. INTRODUCTION 

The T2K (Tokai-to-Kamioka) experiment is a long 
baseline neutrino oscillation experiment [T] whose pri¬ 
mary goal is a precise measurement of the neutrino os¬ 
cillation parameters via the appearance of electron neu¬ 
trinos and the disappearance of muon neutrinos [5]. An 
almost pure intense muon neutrino beam is produced at 
J-PARC (Japan Proton Accelerator Research Complex) 
in Tokai. The proton beam impinges on a graphite tar¬ 
get to produce charged pions, which are focused by three 
magnetic horns [3]. The pions decay mainly into muon- 
muon-neutrino pairs during their passage through the 96- 
meter decay volume. The neutrinos are measured by the 
near detectors (INGRID [3] and ND280 [5HS]) on the J- 
PARC site and the far detector (Super-Kamiokande [TO]) 
in Kamioka, located 295 km away from J-PARC. 

A precise neutrino oscillation measurement requires 
good knowledge of neutrino interaction cross-sections. 
The neutrino charged current quasi-elastic (CCQE) scat¬ 
tering is especially important for T2K because it is used 
as the signal mode for the T2K neutrino oscillation mea¬ 
surement. The I'fj, CCQE cross-section on carbon was 
measured by MiniBooNE [TT], SciBooNE [T2], NOMAD 
[I3], MINERz/A [14] and LSND [15]. The Llewellyn 
Smith formalism m using the relativistic Fermi gas 
model m is generally used to describe CCQE scattering 
with a neutron in the nucleus. However, this approach 
does not provide a good description of existing data. Sev¬ 
eral modifications to the model have been proposed to 
account for the discrepancies, but none of them has yet 
achieved general acceptance. One of the more promising 
approaches involves the introduction of neutrino inter¬ 
actions with two or more nucleons via so-called meson 
exchange current (MEC) into the neutrino interaction 
model [I8ll24j . On the other hand, uncertainties in the 
nuclear model are also regarded as a possible cause of 
the discrepancy [5S]. In order to resolve the puzzle, ad¬ 
ditional CCQE cross-section measurements are required. 

In this paper, we present a measurement of the 
CCQE cross-section on carbon at neutrino energies 
around 1 GeV using the INGRID detector. The CCQE 
signal is defined as the conventional two-body interaction 
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with a single nucleon. We selected the CCQE candidate 
events in INGRID and estimated the CCQE cross-section 
by subtracting background and correcting for selection 
efficiency based on the NEUT neutrino interaction gen¬ 
erator [26]. The CCQE cross-section was estimated as¬ 
suming two different nuclear models, and with and with¬ 
out the multi-nucleon interaction in order to check their 
effects on the cross-section result. 

T2K collected data corresponding to 7.32 x 10^° pro¬ 
tons on target (POT) during the hve run periods listed 
in Table]!] For this cross-section measurement, data from 
Run 2, 3c and 4 are used. The total data set for the cross- 
section measurement corresponds to 6.04 x 10^° POT. 

TABLE I: T2K data-taking periods and integrated 
protons on target (POT). Data of Run 1, 3b, 5a, 5b 
were not used for the cross-section measurement. 


Run Dates 

period 

(Run 1) Jan. 2010 — Jun. 2010 
Run 2 Nov. 2010 — Mar. 2011 
(Run 3b) Mar. 2012 

Run 3c Apr. 2012 — Jun. 2012 
Run 4 Oct. 2012 - May. 2013 
(Run 5a) May. 2014 — Jun. 2014 
(Run 5b) Jun. 2014 


Horn Integrated 
current POT 
250kA 0.32 x lO""" 
250kA 1.11 X 10^° 
205kA 0.21 X 10^“ 
250kA 1.37 X 10^° 
250kA 3.56 x 10^“ 
250kA 0.24 x 10^° 
-250kA 0.51 X 10^° 


The remainder of this paper is organized as follows: 
Details of the INGRID detector and Monte Garlo simula¬ 
tions are explained in Sec. [n]and[ml respectively. Section 


IV summarizes the CCQE event selection. The analysis 


method of the cross-section measurement is described in 
Sec. [V] Section jVIj describes the systematic errors. The 


results and conclusions are given in Sec. VII and VIII 
respectively. 


II. DETECTOR CONFIGURATION 

The INGRID (Interactive Neutrino GRID) detector is 
an on-axis neutrino detector located 280 m downstream 
of the proton target while ND280 and Super-Kamiokande 
are located 2.5° off the beamline axis. It consists of 16 
identical standard modules and an extra module called 
the Proton Module. 

The main purpose of the standard modules is to mon¬ 
itor the neutrino beam direction. Each of the modules 
consists of nine iron target plates and eleven tracking 
scintillator planes. 

In contrast, the Proton Module was developed specifi¬ 
cally for the neutrino cross-section measurement. It is a 
fully-active tracking detector which consists of 36 track- 
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ing layers surrounded by veto planes to reject charged 
particles coming from outside of the modules. The track¬ 
ing layers also serve as the neutrino interaction tar¬ 
get. The total target mass in the fiducial volume is 
303 kg. Seven of the 16 standard modules are horizon¬ 
tally aligned, and the Proton Module is placed in front 
of the central module of them. In this cross-section mea¬ 
surement, the Proton Module is used as the neutrino in¬ 
teraction target, and the seven standard modules located 
downstream of the Proton Module are used as the muon 
detector. The schematic view of the Proton Module and 
the standard modules, and an event display of an MC 
CCQE event are shown in Figs. and respectively. 

The standard module and the Proton Module use a 
common readout system. Scintillation light is collected 
and transported to a photodetector with a wavelength 
shifting fiber (WLS fiber) which is inserted in a hole at 
the center of the scintillator strip. The light is read out by 
a Multi-Pixel Photon Counter (MPPC) EH HI] attached 
to one end of the WLS fiber. The integrated charge and 
timing information from each MPPC is digitized by the 
Trip-t front-end board (TFB) [21]. The integration cycle 
is synchronized with the neutrino beam pulse structure. 
Details of the components and the basic performance of 
the INGRID detector are described in Ref. |3]. 



Proton Module 
Standard modules 


I 


Proton Module Standard module 



FIG. 2: Event display of an MC CCQE event in the 
Proton Module. 


state particles taking into account the intranuclear inter¬ 
actions of hadrons; and a detector response simulation 
which reproduces the final-state particles’ motion and in¬ 
teraction with material, the scintillator light yield, and 
the response of the WLS fibers, MPPCs, and front-end 
electronics. 


A. Neutrino beam prediction 

To predict the neutrino fluxes and energy spec¬ 
tra, a neutrino beam Monte Carlo simulation, called 
JNUBEAM |30|, was developed based on the GEANT3 
framework m- We compute the neutrino beam 
fluxes starting from models (FLUKA2008 [321 1331 ^'^d 
GCALOR [33]) and tune them using existing hadron pro¬ 
duction data (NA61/SHINE [35l]36], Eichten et al. [37] 
and Allaby et al. [55]). Since we use only the beam data 
with the 250kA horn current, the horn current in the sim¬ 
ulation is fixed at 250kA. The predicted neutrino energy 
spectra at the center of INGRID are shown in Fig. [^ 
Energy spectra 10 m upstream of INGRID are predicted 
with the same procedure in order to simulate the back¬ 
ground events from neutrino interactions in the walls of 
the experimental hall. 


FIG. 1: Schematic view of the Proton Module and the 
standard modules. 


III. MONTE CARLO SIMULATION 

The INGRID Monte Carlo (MC) simulation consists 
of three main parts: a simulation of the neutrino beam 
production, which predicts the neutrino flux and energy 
spectrum of each neutrino flavor; a neutrino interaction 
simulation, which is used to calculate the neutrino in¬ 
teraction cross-sections and the kinematics of the final 


B. Neutrino interaction simulation 

Neutrino interactions with nuclear targets are simu¬ 
lated with the NEUT program library [26] . Both the pri¬ 
mary neutrino interactions in nuclei and the secondary 
interactions of the hadrons in the nuclear medium are 
simulated in NEUT. NEUT uses the Llewellyn Smith 
formalism [16] with the relativistic Fermi gas model for 
quasi-elastic scattering, the Rein-Sehgal model [521130] 
for resonant meson production and coherent tt produc¬ 
tion and GRV98 (Gluck-Reya-Vogt-1998) [ 31 ] parton dis¬ 
tributions with Bodek-Yang modifications [321 [ 33 ] for 
deep inelastic scattering. For the measurement presented 
in this paper, the final state interactions of the nude- 
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FIG. 3: Neutrino energy spectrum for each neutrino 
species at the Proton Module with the 250kA horn 
current predicted by JNUBEAM. 


ons in nuclei are important, because we use information 
about protons produced in the interaction. In NEUT, 
both elastic scattering and pion production are consid¬ 
ered. The differential cross-sections were obtained from 
nucleon-nucleon scattering experiments |44] . For pion 
production, the isobaric nucleon model [15] is used. The 
TT-less A decay, which is the interaction of the A within 
the target nucleus prior to decay into a pion and a nu¬ 
cleon, is also an important process because the reso¬ 
nant interaction event is often misidentified as the CCQE 
event due to this interaction. We use the data from 
Refs. [461147] to determine the probability of the num¬ 
ber of emitted nucleons and the kinematics of the two- 
nucleon emission in the 7r-less A decay. For emission 
of three or more nucleons, the nucleons are isotropically 
emitted. The simulations in NEUT are described in more 
detail in Refs. [5^ US]. Figure shows the neutrino- 
nucleus cross-sections per nucleon divided by the neu¬ 
trino energy predicted by NEUT. Additionally, a CCQE 
cross-section prediction by a different neutrino interac¬ 
tion simulation package, GENIE [13], is used for com¬ 
parison. GENIE also uses the Llewellyn Smith formal¬ 
ism with the relativistic Fermi gas model. However, the 
nominal value of the axial mass [50] differs from that 
in NEUT (1.21 GeV/c^ for NEUT and 0.99 GeV/c^ for 
GENIE). In addition, GENIE incorporates short range 
nucleon-nucleon correlations in the relativistic Eermi gas 
model and handles kinematics for off-shell scattering ac¬ 
cording to the model of Bodek and Ritchie m while 
NEUT uses the Smith-Moniz model [T7] . 


C. INGRID detector response simulation 

The INGRID detector simulation was developed using 
the Geant4 framework [52]. It models the real detector 
structures (geometries, materials). The structure of the 
walls of the experimental hall is also modeled to sim¬ 
ulate background events from neutrino interactions in 



FIG. 4: Neutrino-nucleus cross-sections per nucleon of 
carbon nucleus divided by the neutrino energy, as 
predicted by NEUT. 


the walls. The particles’ motion and interactions with 
the materials are simulated, and the energy deposit of 
each particle inside the scintillator is stored. Simulations 
of hadronic interactions are performed with the QGSP 
BERT physics list [53]. The energy deposit is converted 
into a number of photons. Quenching effects of the scin¬ 
tillation are modeled based on Birks’ law [Ml [55]. The 
effect of collection and attenuation of the light in the scin¬ 
tillator and the WLS fiber is modeled based on the results 
of electron beam irradiation tests. The non-linearity of 
the MPPG response is also taken into account, since the 
number of detectable photoelectrons is limited by the 
number of MPPG pixels. The number of photoelectrons 
is smeared according to statistical fluctuations and elec¬ 
trical noise. The dark count of the MPPCs is added with 
a probability calculated from the measured dark rate. 
Because the response of the ADGs on front-end electron¬ 
ics is not linear, its response is modeled based on the 
results of a charge injection test. 

IV. EVENT SELECTION 
A. CC event selection 

As the first step, GG interaction events in the Proton 
Module are reconstructed and selected as follows where 
the coordinates shown in Fig. [^are used. 

1. Two-dimensional XZ and YZ tracks in the Pro¬ 
ton Module and the standard modules are re¬ 
constructed independently from hit information in 
each plane using a specially developed reconstruc¬ 
tion algorithm. 

2. Track matching between the Proton Module tracks 
and the standard module tracks are performed. 

3. Three-dimensional tracks are searched for among 
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pairs of two-dimensional XZ tracks and YZ tracks 
by matching the Z positions of the track edges. 

4. The neutrino interaction vertices are searched for 
by looking at the upstream edges of the three- 
dimensional tracks. 

5. Events within ±100 nsec from the expected timing, 
which is calculated from the timing of the pulsed 
primary proton beam, are selected. 

6. Events which have a hit in a veto plane at the up¬ 
stream position extrapolated from a reconstructed 
track are rejected. The events rejected by the front 
veto plane are identified as beam induced muon 
backgrounds. The number of neutrino interactions 
on the walls in the MC simulation is normalized by 
the observed number of the beam induced muon 
backgrounds. 

7. Events having a vertex outside the fiducial volume, 
which is defined as within ±50 cm from the module 
center in the X and Y directions, are rejected. 

Further details of this neutrino event selection is written 
in Ref. [56]. After this selection, CCQE events make up 
37.67% of the MC sample. To increase the selection pu¬ 
rity for CCQE events, additional cuts are applied based 
on the number of reconstructed tracks, the dE/dx par¬ 
ticle identification variable and the reconstructed event 
kinematics. Then the selected CCQE candidate events 
are classified according to the neutrino energy. 


B. Number of tracks 

1. Number of tracks from the vertex 

The CCQE interaction produces two particles inside 
the target nucleus, a muon and a proton. However, the 
proton undergoes final state interactions (ESI) in the 
residual nucleus and does not always escape unaltered. 
In addition, the proton may not be reconstructed in the 
detector due to its short range. Thus, events with either 
one or two reconstructed tracks coming from the vertex 
are selected (Fig. [^. Hereafter, the events with one and 
two reconstructed tracks are referred to as the one-track 
sample and the two-track sample, respectively. The anal¬ 
ysis must account for the above effects to produce the 
final cross-section result from each sample. 


2. Number of matched tracks 

When selecting CC events, events with at least one 
matched track between the Proton Module and the stan¬ 
dard module are accepted in order to select those con¬ 
taining a long muon track. About 7% of the selected 
events in the two-track sample have two matched tracks. 



FIG. 5: The number of reconstructed tracks from vertex 
in the Proton Module. The colored histograms are the 
MC predictions divided by the neutrino type and 
interaction type. The neutrino interaction events in the 
Proton Module are absolutely normalized to POT, and 
the background events from outside are normalized to 
beam induced muon backgrounds. 


Simulations indicate that the second matched track is 
usually a pion from a CC-nonQE interaction such as CC 
resonant pion production, CC coherent pion production 
etc. Thus, events with exactly one matched track are se¬ 
lected for the two-track sample (Fig. |^. Hereafter, the 
matched track is referred to as the first track and the 
remaining track as the second track. 



FIG. 6: The number of matched tracks between the 
Proton Module and the standard module for the two 
track sample. 
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C. Particle identification 

1. Definition of the muon confidence level (MuCL) variable 

Particle identification (PID) based on dE/dx informa¬ 
tion is applied on both the one-track sample and the 
two-track sample. The dE/dx for each scintillator plane 
is calculated from the light yield divided by the path 
length of the track in the scintillator where the light at¬ 
tenuation in the WLS fiber is corrected. The first step 
of the particle identification is to estimate a confidence 
level that a particle is a muon on a plane-by-plane ba¬ 
sis. The confidence level at each plane is defined as the 
fraction of events in the expected dE/dx distribution of 
muons above the observed dE/dx value because muons 
have a lower dE/dx value than protons. The expected 
dE/dx distribution of muons is obtained from the beam- 
induced muon backgrounds which are mainly created by 
the neutrino interactions in the walls of the detector hall. 
The cumulative distribution function of the muon dE/dx 
distribution corresponds to the confidence level. The cal¬ 
culated confidence level at the f-th plane as a function of 
dE/dx is referred to as CLi. 

The next step is to combine the confidence levels (CLi) 
obtained from all the planes penetrated by the track to 
form a total confidence level. In the case where the 
track penetrates only two planes, the procedure to com¬ 
bine the two confidence levels, CLi and CL 2 , is as fol¬ 
lows. Assuming the confidence levels to be independent 
of one another, the combined probability is the product, 
P = CLi X CL 2 . In the a:?/-plane of the two confidence 
levels, the hyperbola xy = P gives such a combined prob¬ 
ability and the unified muon confidence level, MuCL, is 
the fraction of possible x,y values that give xy < P. 
MuCL is expressed as 



= P(I-lnP). (1) 

In analogy with the two-plane case, the muon confidence 
level combined from n planes is expressed as 

MuCL = P X ^ ^ ^ , P = n (2) 

i—0 z—1 

When a Proton Module track is matched with a standard 
module track, the standard module track is also used to 
make the MuCL. 


2. Separation of muon-like tracks and proton-like tracks 

The last step of the particle identification is to distin¬ 
guish the tracks using the MuCL. In this analysis, tracks 
whose MuCL are more than 0.6 are identified as muon¬ 
like and those less than 0.6 are identified as proton-like. 


The probability of misidentifying a muon track (a pro¬ 
ton track) as proton-like (muon-like) in the MC simula¬ 
tion is 12.5% (10.9%). Most pion tracks are identified 
as muon-like, since the mass of the pion is similar to the 
muon mass. For the one-track sample, events having a 
muon-like track (Fig. are selected as the CCQE en¬ 
hanced sample. For the two-track sample, the first track 
is required to be muon-like and the second track to be 
proton-like (Fig. |^|^. 




FIG. 8: MuCL distributions of the first track for the 
two-track sample. 


D. Kinematic cut 

In addition, two kinematic cuts are applied to the two- 
track sample. These cuts use two angles called the copla¬ 
narity angle and the opening angle, defined as shown in 

Fig.[^ 



































FIG. 9: MuCL distributions of the second track for the 
two-track sample. 


FIG. II: Goplanarity angle distribution following the 
PID cut for the two track sample. 


2. Opening angle cut 



FIG. 10: Definition of the coplanarity angle and the 
opening angle. 


1. Coplanarity angle cut 


Since CCQE events are (quasi) two-body scattering in¬ 
teractions, all the tracks in a GGQE event (an incident 
neutrino track, a scattered muon track and a scattered 
proton track) are expected to lie in the same plane if the 
effects of the proton re-scatterings and the Fermi mo¬ 
mentum of the target nucleons are neglected. To quan¬ 
tify this, the coplanarity angle is defined as the angle 
between the two reconstructed three-dimensional tracks 
projected to the XY plane, where the XY plane is per¬ 
pendicular to the neutrino beam axis (Fig. 10). When the 
three tracks are precisely coplanar, the coplanarity angle 
is 180°. Thus, events with a coplanarity angle above 150° 
are selected (Fig. 111. 


The opening angle is defined as the angle between 
the two reconstructed three-dimensional tracks (Fig. 101. 
The opening angle tends to be large in CCQE interac¬ 
tions, because in the center of mass frame the two final 
particles are produced back to back. Thus, events with 
an opening angle above 60° are selected (Fig. [^. The 
results of the event selection so far are summarized in 
Table m 



FIG. 12: Opening angle distribution following the PID 
and coplanarity angle cuts for the two track sample. 


E. Energy classification 


We aim to measure the CCQE cross-section in the 
low energy region (^1 GeV) and the high energy region 
(^2 GeV) separately. Thus, an energy classification is 
applied to the CCQE enhanced samples to select subsam¬ 
ples enriched in high-energy and low-energy events re¬ 


spectively. The classification criterion is shown in Fig. 13 
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TABLE II: The number of events passing each CCQE 
selection step. The efficiency is defined as the number 
of selected CCQE events divided by the number of 
CCQE interactions in the FV. The purity is defined 
as the ratio of the selected CCQE events to the total 
selected events. 


Selection 

Data MC Efficiency Purity 

One track from vertex 

12896 1.23x10^ 

25.5% 

60.1% 

Particle identification 

9059 8.75x10^ 

23.0% 

76.1% 

Two tracks from vertex 

14479 1.47x10"* 

20.1% 

39.4% 

One matched track 

13436 1.37x10^ 

19.6% 

41.3% 

Particle identification 

7981 8.32x10^ 

15.8% 

54.9% 

Kinematic cut 

3832 4.23x10^ 

12.2% 

83.5% 


Events with a muon candidate track which penetrates all 
the standard module iron layers are selected as the high- 
energy sample, while events with a muon candidate track 
which stops in the standard module are selected as the 
low-energy sample. Other events with a muon candidate 
track which escapes from the side of the standard module 
are not used in this analysis. Figures [T^ andshow the 
neutrino energy spectra of the CCQE enhanced samples 
before and after applying the energy classification. Most 
of the CCQE events at neutrino energies below I.O GeV 
(above 1.5 GeV) are rejected by the high energy selection 
(the low energy selection). Figure 16 shows the selec¬ 


tion efficiency of each sample as a function of the muon 
momentum and angle. Each sample covers the different 
muon kinematic region. 


Proton Module 


Standard module 



FIG. 13: Event display of penetrating, stopped and 
side-escaped events. 


F. Event pileup correction 

The T2K neutrino beam is pulsed. Each pulse has 
an eight-bunch structure, and each bunch has a width 
of 58 ns. When a track from a neutrino event piles up 
with a track from another neutrino event in the same 



FIG. 14: True neutrino energy spectra of the GGQE 
enhanced sample before and after applying the energy 
classification for the one-track sample in the MG 
simulation. The spectrum before the energy 
classification includes the side-escaped events. 



FIG. 15: True neutrino energy spectra of the GGQE 
enhanced sample before and after applying the energy 
classification for the two-track sample in the MG 
simulation. The spectrum before the energy 
classification includes the side-escaped events. 


beam bunch, vertices may fail to be reconstructed. Be¬ 
cause this results in the loss of events, this event-pileup 
effect needs to be corrected for. The event-pileup effect 
is proportional to the beam intensity. Hence, the correc¬ 
tion factor is estimated as a linear function of the beam 
intensity, where the slope of the linear function is esti¬ 
mated from beam data as follows. First, the beam data 
is categorized into subsamples according to the beam in¬ 
tensity. In each subsample, all hits in INGRID from two 
beam bunches are summed together to make one new 
pseudo beam bunch. This procedure effectively doubles 
the beam intensity observed by INGRID. A slope is esti¬ 
mated from the number of selected events in an original 
beam bunch and a pseudo beam bunch for each subsam¬ 
ple. The slopes estimated from all subsamples are consis¬ 
tent with each other, and the average value of this slope 
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Muon momentum (GeV/c) 


(a) High energy one-track sample 



Muon momentum (GeV/c) 


(c) High energy two-track sample 
FIG. 16: Selection efficiency of each sample as 



Muon momentum (GeV/e) 


(b) Low energy one-track sample 



o t.-I -Iq 

0 0.5 1 1.5 2 2.5 3 3.5 


Muon momentum (GeV/c) 
(d) Low energy two-track sample 

a function of the muon momentum and angle. 


is used for the correction. This event pileup correction 
is applied bunch-by-bunch using the slope and POT per 
bunch which corresponds to the relevant beam intensity. 
The event pileup correction gives 0.3-0.7% difference in 
the number of selected event in each sample. 


V. ANALYSIS METHOD 


We estimate the CCQE cross-sections in high and low 
energy regions, which are defined as above 1.0 GeV and 
below 1.5 GeV using the high and low energy samples. 
The average energies of the neutrino flux in the high and 
low energy regions are 1.94 GeV and 0.93 GeV, respec¬ 
tively. The GGQE cross-section is calculated from the 
number of selected CCQE candidate events by subtract¬ 
ing background and correcting for selection efficiency: 


CTCCQE = 


Nael — NbG 


(3) 


where Wei is the number of selected CCQE candidate 
events from real data, Nbg is the number of selected 


background events predicted by the MC simulation, </> is 
the integrated flux, T is the number of target neu¬ 
trons, and e is the detection efficiency for CCQE events 
predicted by the MC simulation. The flirx; (j) is inte¬ 
grated in each energy region, and £ is calculated for 
CCQE events in each energy region. CCQE events as¬ 
signed to the wrong subsample, i.e. high-energy events 
in the low-energy subsample and vice versa, are regarded 
as background. The background events for this analy¬ 
sis consist of CC-nonQE events, NC events, W events, 
z/e events, external background events created by neu¬ 
trino interactions in the material surrounding the de¬ 
tector, and CCQE events assigned to the wrong energy 
sample. Furthermore, the CCQE cross-section in each 
energy region is estimated from the one-track sample, 
two-track sample, and combined sample, separately. The 
one-track sample (the two-track sample) has an enhanced 
content of low (high) energy protons from CCQE interac¬ 
tions as shown in Fig. |17| Therefore, we can cross-check 
the CCQE cross-section results from the different phase 
spaces. 





















11 



FIG. 17: Distributions of the true kinetic energy of 
protons from the CCQE events in the one-track sample 
and the two-track sample in the MC simulation. The 
cutoff around 0.02 GeV stems from the effect of Pauli 
blocking. 


VI. SYSTEMATIC ERRORS 

Uncertainties on A^bGj T, and e are sources of sys¬ 
tematic errors on the cross-section results. The sources 
of systematic error can be categorized into three groups: 
those from the neutrino flux prediction, the neutrino in¬ 
teraction model including intra-nuclear interactions, and 
the detector response. 


A. Neutrino flux uncertainties 

The neutrino flux uncertainty sources can be sepa¬ 
rated into two categories: hadron production uncertain¬ 
ties and T2K beamline uncertainties. The uncertain¬ 
ties on hadron production are mainly driven by the 
NA61/SHINE measurements [351I3S] and the Eichten and 
Allaby data |37l[38], and constitute the dominant compo¬ 
nent of the flux uncertainty. They include the uncertain¬ 
ties on the production cross-section, the secondary nu¬ 
cleon production, the pion production multiplicity, and 
the kaon production multiplicity. The second category 
of flux uncertainties is associated with inherent uncer¬ 
tainties and operational variations in the beamline con¬ 
ditions, including uncertainties in the proton beam posi¬ 
tion, the beam-axis direction, the absolute horn current, 
the horn angular alignment, the horn field asymmetry, 
the target alignment, and the proton beam intensity. The 
method of estimating these flux uncertainties is described 
in Ref. [30]. To evaluate the systematic error from the 
flux uncertainties, the flux is varied using a covariance 
matrix based on the flux uncertainty. This is repeated 
for many toy data sets, and the ±lcr of the change in 
the cross-section result is taken as the systematic error 
associated with the neutrino flux. The systematic error 
is 11-17%, which is the dominant error in this measure¬ 


ment. 


B. Neutrino interaction uncertainties 

We use a data-driven method to calculate the neutrino 
interaction uncertainties, where the NEUT predictions 
are compared to available external neutrino-nucleus data 
in the energy region relevant for T2K. We fit some param¬ 
eters of the models implemented in NEUT, and introduce 
ad hoc parameters, often with large uncertainties, to take 
into account remaining discrepancies between NEUT and 
the external data [TTJ [TS] jM] [57H65] . 

The model parameters include axial mass values for 
quasi-elastic scattering and meson production via baryon 
resonances, the Fermi momentum, the binding energy, a 
spectral function parameter, and a 7r-less A decay pa¬ 
rameter. The spectral function parameter is introduced 
to take into account the difference between the rela¬ 
tivistic Fermi gas nuclear model, which is the standard 
NEUT model, and the more sophisticated spectral func¬ 
tion model |66j . which is expected from electron scat¬ 
tering data [SZIISHI to be a better representation of the 
nuclear model. The implemented ad hoc parameters in¬ 
clude cross-section normalizations. In addition, uncer¬ 
tainties on the final state interactions of the pions and 
nucleons with the nuclear medium are included. Ta¬ 
ble mil shows the nominal values and uncertainties on 
these parameters. Systematic errors from the nuclear 
model (spectral function, Fermi momentum and binging 
energy) were found to be comparatively large because the 
CCQE interaction in the few GeV region is sensitive to 
the nuclear model. Further details about these uncer¬ 
tainties are described in Refs UHl IM| • Systematic errors 
due to these parameters are estimated from variations of 
the cross-section results when these parameters are var¬ 
ied within their uncertainties. As a cross check, we also 
estimated the selection efficiencies using the GENIE neu¬ 
trino interaction generator and confirmed that they were 
consistent with those using NEUT when the axial vector 
mass is set to the same value as GENIE, 0.99 GeV. 


C. Detector response uncertainties 

The uncertainty of the target mass measurement, 
0.25%, is taken as the systematic error on the target 
mass. The variation in the measured MPPC dark rate 
during data acquisition, 11.52 hits/cycle, where ‘cycle’ 
denotes the integration cycle synchronized with the neu¬ 
trino beam pulse structure, is taken as the uncertainty on 
the MPPG dark rate. The discrepancy between the hit 
detection efficiency measured with beam-induced muon 
backgrounds and that of the MC simulation, 0.21%, is as¬ 
signed as the uncertainty in the hit detection efficiency. 
The uncertainty of the light yield is evaluated by us¬ 
ing the beam-induced muon backgrounds as a control 
sample; in addition, the uncertainty of the scintillator 
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TABLE III: The nominal values of and the uncertainties 
on the neutrino interaction model parameters. The 
first, second, and third groups represent the model 
parameters, the ad hoc parameters, and the final state 
interaction parameters respectively. The 7r-less A decay 
parameter and the final state interaction parameters 
vary the probabilities of these interactions. Nominal 
values of 0 or 1 mean that the effect or the 
normalization is not implemented or is implemented by 
default, respectively. 


tween data and MC for the subsamples is taken as the 
systematic error. The systematic errors from all event 
selections are evaluated by varying each selection thresh¬ 
old. The maximum difference between the data and MC 
for each selection threshold is taken as the value of each 
systematic error. Among the systematic errors from the 
detector response, the largest contributions (about 2% 
each) are those from the light yield and the secondary 
interactions. 


Parameter 

Nominal 

Error 


1.21GeV 

16.53% 


1.21GeV 

16.53% 

TT-less A decay 

0.2 

20% 

Spectral function 

O(off) 

100% 

Fermi momentum 

217MeV/c 

13.83% 

Binding energy 

25MeV 

36% 

CCItt norm. {E^, <2.5GeV) 

1 

21% 

CCItt norm. (E„ >2.5GeV) 

1 

21% 

GC coherent tt norm. 

1 

100% 

GC other shape 

O(off) 

40% at 1 GeV 

NCItt® norm. 

1 

31% 

NG coherent tt norm. 

1 

30% 

NCItt^ norm. 

1 

30% 

NC other norm. 

1 

30% 

TT absorption 

1 

50% 

TT charge exchange (low energy) 

1 

50% 

TT charge exchange (high energy) 

1 

30% 

TT QE scattering (low energy) 

1 

50% 

TT QE scattering (high energy) 

1 

30% 

TT inelastic scattering 

1 

50% 

Nucleon elastic scattering 

1 

10% 

Nucleon single tt production 

1 

10% 

Nucleon two tt production 

1 

10% 


quenching is taken into account based on the uncertainty 
of the Birks’ constant (0.0208±0.0023 cm/MeV). The re¬ 
lations between these quantities and the cross-section re¬ 
sults are estimated by MC simulation, and the resulting 
variations on the calculated cross-sections are assigned as 
systematic errors. The event pileup correction factor has 
uncertainties which come from the statistics of the beam 
data and the MPPC dark count in the estimation of the 
correction factor. The systematic error from these uncer¬ 
tainties is estimated assuming the highest beam intensity 
achieved in beam operation so far. There is about a 
35% discrepancy between the beam-induced muon back¬ 
ground rate estimated by the MC simulation and that 
measured from the data. The change in the background 
contamination fraction from this discrepancy is taken as 
the systematic error for the beam-related background. 
The cosmic-ray background was found to be very small 
from the out-of-beam timing data. The systematic er¬ 
ror on the track reconstruction efficiency is estimated by 
comparing the efficiency for several subsamples between 
the data and the MC simulation. The standard deviation 
of the difference of the track reconstruction efficiency be¬ 


D. Summary of the systematic errors 


The total systematic errors, calculated from the 
quadrature sum of all the systematic errors, on the 
CCQE cross-section from the one-track sample, the two- 


track sample and the combined sample for the high en- 

prcv rpcinn arp +15-95% +16.97% +15.06% ,.pqT-,pp_ 

region are _X 2 . 97 %’ — 14 . 04 % — 12 . 44 %’ respec 

tively, and those for the low energy region are 

summarizes 


-18.86% 


and _45 49 %, respectively. Table 


IV 


the breakdown of the systematic errors on the CCQE 
cross-section measurement from the combined sample. 


The neutrino flux error is the dominant systematic er¬ 


ror. 


VII. RESULTS 


A. Results with the T2K default interaction model 


The results of the CCQE cross-section measurement 
from the different subsamples are summarized in Table 
[^(a). The measured CCQE cross-sections from the com¬ 
bined sample are 

(7ccQE(1.94GeV) = (11.95 ± 0.19(stat.)ti;®°(si/st.)) 

X 10“^® cm^/neutron (4) 

^TccQE(0.93GeV) = (10.64 ± 0.37(stat.)l?;°3(sj/st.)) 

X 10“^®cm^/neutron, (5) 


at mean neutrino energies of 1.94 GeV and 0.93 GeV (for 
the high and low energy regions), respectively. We quote 
these values as our primary result. The NEUT and GE¬ 
NIE predictions of the CCQE cross-sections on carbon for 


the high and low energy regions are shown in Table VI 


The difference in the predictions from NEUT and GENIE 
is attributable primarily to the difference in the nominal 
value. The results of the measurements are consis¬ 
tent within 2(7 with both predictions; however, in the low- 
energy region the cross-section results from the one-track 
and two-track samples differ by just under 2cr, as shown 
in Table [Vnl The cross-section results are shown in Fig. 
[TS] together with the predictions and the measurements 
of other experiments, and the reconstructed muon angle 
distribution of each CCQE candidate sample is shown in 
Fig. The distributions of other kinematic variables 
are summarized in the supplemental material CO]. 
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TABLE IV: Summary of the systematic errors on the CCQE cross-section measurement from the combined sample. 

Negative and positive values represent — Icr and -|-lcr errors. 


Item 

High energy region 

Low energy region 

Neutrino flux 

-11.01% -b 13.61% 

-13.57% ± 17.04% 


-0.89% -b 2.25% 

-0.08% ± 0.39% 

MREs 

-0.92% -b 1.31% 

-0.82% ± 1.10% 

CCItt normalization {Ei, <2.5 GeV) 

-0.55% -b 0.50% 

-3.71% ± 3.59% 

CCItt normalization {Ei, >2.5 GeV) 

-2.69% -b 2.69% 

-1.88% ±1.83% 

GC coherent tt normalization 

-1.40%-b 1.38% 

-1.73% ± 1.71% 

GC other E^ shape 

-0.86% -b 0.85% 

-0.11% ±0.09% 

NGItt® normalization 

-0.65% -b 0.65% 

-0.40% ± 0.40% 

NG coherent tt normalization 

-0.10%-b 0.10% 

-0.09% ± 0.09% 

NGItt^ normalization 

-0.47% -b 0.47% 

-0.46% ± 0.45% 

NG other normalization 

-0.33% -b 0.31% 

-0.75% ± 0.74% 

TT-less A decay 

-0.54%-b 2.10% 

-1.60% ± 3.34% 

Spectral function 

-2.01% -b 0.00% 

-0.00% ± 1.21% 

Fermi momentum 

-1.67%-b 2.22% 

-3.71% ± 4.43% 

Binding energy 

-0.44% -b 0.65% 

-1.24% ±1.42% 

Pion absorption 

-0.20% -b 0.81% 

-0.80% ±1.20% 

Pion charge exchange (low energy) 

-0.15%-b 0.18% 

-0.22% ± 0.28% 

Pion charge exchange (high energy) 

-0.11%-b 0.13% 

-0.11% ±0.11% 

Pion QE scattering (low energy) 

-0.66% -b 0.71% 

-0.84% ± 0.79% 

Pion QE scattering (high energy) 

-0.04% -b 0.03% 

-0.09% ± 0.09% 

Pion inelastic scattering 

-0.05% -b 0.04% 

-0.29% ± 0.25% 

Nucleon elastic scattering 

-0.25% -b 0.21% 

-0.29% ± 0.21% 

Nucleon single tt production 

-0.15%-b 0.11% 

-0.60% ± 0.51% 

Nucleon two tt production 

-0.57% -b 0.42% 

-0.01% ± 0.01% 

Target mass 

±0.31% 

±0.38% 

MPPC dark noise 

±0.03% 

±0.08% 

Hit efflciency 

±0.84% 

±0.41% 

Light yield 

±1.47% 

±2.22% 

Event pileup 

±0.02% 

±0.06% 

Beam-induced external background 

±0.08% 

±0.35% 

Gosmic-ray background 

±0.00% 

±0.01% 

2D track reconstruction 

±0.67% 

±0.81% 

Track matching 

±0.45% 

±1.13% 

3D tracking 

±0.21% 

±0.15% 

Vertexing 

±0.30% 

±0.43% 

Timing cut 

±0.00% 

±0.00% 

Veto cut 

±0.82% 

±0.64% 

Fiducial volume cut 

±1.55% 

±0.84% 

Secondary interaction 

±2.45% 

±2.37% 

Total 

-12.44% ±15.06% 

-15.49% ± 19.04% 


B. Results with the spectral function model 

The T2K default interaction model uses the relativis¬ 
tic Fermi gas model as the nuclear model. As dis¬ 
cussed above, the spectral function model is more so¬ 
phisticated and is expected to provide a better descrip¬ 
tion of neutrino-nucleus interactions. When O. Benhar’s 
spectral function m is used in the MC simulation for 
the efficiency correction instead of the relativistic Fermi 
gas model, the CCQE cross-section results are slightly 
changed, as shown in Table [V| (b), because of the differ¬ 
ences in the kinematics of the final state particles which 
arise from the differences in the initial nucleon momen¬ 
tum distribution. The cross-section results derived using 


the spectral function model are shown in Fig. [^together 
with the model predictions. The cross-section result in 
the low energy region from the one-track sample (the 
two-track sample) with the spectral function is 5% lower 
(10% higher) than that with the relativistic Fermi gas 
model. As a result, the difference in the cross-section re¬ 
sults between one-track and two-track samples in the low 
energy region becomes smaller, as shown in Table [VTl| in 
which correlations of the systematic errors between the 
cross-section results are correctly treated. This change 
mainly comes from the differences in the final state pro¬ 
ton kinematics between the relativistic Fermi gas model 
and the spectral function, which cause event migrations 
between the one-track sample and the two-track sample. 
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TABLE V: The CCQE cross-sections measured from each sample (xlO ^®cm^). The mean neutrino energies of the 
high and low energy regions are 1.94 GeV and 0.93 GeV, respectively. 


(a) Results with the relativistic Fermi gas model (T2K default interaction model). 


Used sample High energy region 

Low energy region 

One-track sample 12.29 ± 0.22(stat.)i((^'gg(s2/st.) 
Two-track sample 10.98 ± 0.35(stat.)l)) g 4 (s 2 /st.) 
Combined sample 11.95 ± 0.19(stat.)l4 |g(syst.) 

11.63 ± 0.45(stat.)+f gg(syst.) 
8.01 ± 0.64(stat.)l} j4(syst.) 

10.64 ± 0.37{stat.)+l °l{sySt.) 

(b) Results with the spectral function model. 

Used sample High energy region 

Low energy region 

One-track sample 12.46 ± 0.22(stat.)i(4'g2(s2/st.) 
Two-track sample 11.43 ± 0.36(stat.)l4 gg(s 2 /st.) 
Combined sample 12.19 ± 0.19(stat.)l4 gg(syst.) 

11.04 ± 0.43{stat.)1:i%l{syst.) 

8.84 ± 0.70(staL)t};?o(syst.) 
10.51 ± 0.37{stat.)tl °°{^ySt.) 

(c) Results on assuming the multi-nucleon interactions. 

Used sample High energy region 

Low energy region 

One-track sample 10.79 ± 0.22{stat.)'^f^^{syst.) 
Two-track sample 10.28 ± 0.35(stat.)l4 g|(s?/st.) 
Combined sample 10.66 ± 0.19(stat.)l4 g|(syst.) 

10.11 ± 0.45(stat.)i^ "(syst.) 

7.14 ± 0.6A{stat.)tl-lllsyst.) 
9.30 ± 0.37{stat.)+l-ll{syst.) 


TABLE VI: The NEUT and GENIE predictions of the 
flux-averaged GGQE cross-sections on carbon for the 
high and low energy regions. 


High energy region Low energy region 
NEUT ITMxICp^^Vn? 10.34x 10"^“cm^ 

GENIE 9.46xl0"^®cm^ 8.49xl0"®®cm^ 


xlO“ 


a 

o 

o 
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■ NEUT(Mg^=l-21GeV/c^l 

■ GENIE {MSM.99GeV/c^)- 

■ NEUT flux average 
GENIE flux average - 

10 

E, (GeV) 


FIG. 18: The CCQE cross-section results with 
predictions by NEUT and GENIE. Our data point is 
placed at the flux mean energy. The vertical error bar 
represents the total (statistical and systematic) 
uncertainty, and the horizontal bar represent 68% of the 
flux at each side of the mean energy. The SciBooNE, 
MiniBooNE, NOMAD and MINERi/A results are also 
plotted UMM- 


Therefore, this result may indicate that the spectral func¬ 
tion is a better representation of the nuclear model than 
the relativistic Fermi gas model. 


C. Results on assuming multi-nucleon interactions 


The T2K default interaction model does not assume 
the existence of multi-nucleon interactions via the me¬ 
son exchange current. Multi-nucleon interactions gen¬ 
erally produce a lepton and two or more nucleons, but 
are misidentified as GGQE events when the additional 
nucleons are not reconstructed. This means that the 
CCQE cross-section measurement is expected to be sensi¬ 
tive to the existence of multi-nucleon interactions. There¬ 
fore, we also estimate the CCQE cross-sections assuming 
the existence of multi-nucleon interactions. There are 
many multi-nucleon interaction models |I8H24j : we used 
the model proposed by J. Nieves Eoi mi- When the 
multi-nucleon interaction model was introduced into the 
neutrino interaction model, the CCQE model including 
the value of and the relativistic Fermi gas model 
were not changed in order to isolate the effect of the 
multi-nucleon interaction. In this CCQE cross-section 
analysis, the CCQE signal is defined as the conventional 
two-body interaction with a single nucleon. Therefore, 
the multi-nucleon interaction events are defined to be 
background events and are subtracted from the selected 
events as with CC-nonQE events, NC events, etc. The 
CCQE cross-section results assuming the existence of 
multi-nucleon interactions are summarized in Table El 


(c) and shown in Fig. 21 together with the predictions. 
Although the measured CCQE cross-sections assuming 
the existence of the multi-nucleon interaction are 6-13% 
smaller, they are still compatible with the predictions. 















































Number of events Number of events 
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(a) High energy one-track sample 


(b) Low energy one-track sample 




(c) High energy two-track sample (d) Low energy two-track sample 


FIG. 19: Reconstructed muon angle distribution of each CCQE candidate sample. Explanations about the 
predictions with the spectral function model and the multi-nucleon interaction model will be given in Secs. |VII B| 


and VIIC 
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TABLE VII: Ratio of the CCQE cross-section result from the one-track sample to that from the two-track in the 

low energy region. 


Model in MC simulation Ratio of cross-section results 

Relativistic Fermi gas model (T2K default model) 1.45 ± 0.09{stat.)'^Q'2g{syst.) 

Spectral function model 1.25 ± 0.0S{stat.)'^Q'2g{syst.) 

Relativistic Fermi gas model with multi-nucleon interactions 1.42 ± 0.09{stat.)t.Q'^Ji{syst.) 



FIG. 20: The CCQE cross-section results when the 
spectral function is used as the nuclear model. 



FIG. 21: The CCQE cross-section results on the 
assumption of the existence of the multi-nucleon 
interaction. The relativistic Fermi gas model is used as 
the nuclear model. 


VIII. CONCLUSIONS 

We have reported a CCQE cross-section measurement 
using the T2K on-axis neutrino detector, INGRID. We 
have selected one-track and two-track samples of 
CCQE scattering in the Proton Module. From the num¬ 
ber of selected events, the CCQE cross-sections on car¬ 
bon at mean neutrino energies of 1.94 GeV and 0.93 GeV 
have been measured. The cross-section analysis was per¬ 
formed using three different neutrino interaction models: 
relativistic Fermi gas models with and without multi¬ 
nucleon interactions, and a spectral function model with¬ 
out multi-nucleon interactions. Although these results 
are compatible with the model predictions, we found that 
the CCQE cross-section results are dependent on the nu¬ 
clear model and the existence of multi-nucleon interac¬ 
tions at the 10% level. There is some indication, based 
on consistency between the one-track and two-track sam¬ 
ples, that the event kinematics are better described by 
the spectral function model. The data related to this 
measurement can be found electronically in m- 
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